Anomalous suppression of the shot noise in a nanoelectromechanical system 
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In this paper we report a relaxation-induced suppression of the noise for a single level quantum 
dot coupled to an oscillator with incoherent dynamics in the sequential tunneling regime. It is 
shown that relaxation induces qualitative changes in the transport properties of the dot, depending 
on the strength of the electron-phonon coupling and on the applied voltage. In particular, critical 
thresholds in voltage and relaxation are found such that a suppression below 1/2 of the Fano factor 
is possible. Additionally, the current is either enhanced or suppressed by increasing relaxation, 
depending on bias being greater or smaller than the above threshold. These results exist for any 
strength of the electron-phonon coupling and are confirmed by a four states toy model. 

PACS numbers: 73.50.Td,73.23.-b,85.85.+j 



I. INTRODUCTION 



In the last years, nanoelectromechanical systems 
(NEMS) have been a hot research topic both from the 
theoretical and the experimental point of view. i£ Com- 
bining electronic and mechanical degrees of freedom, 
NEMS have potentially important applications as fast 
and ultra-sensitive detectors; 3 ^^ as well as being in- 
teresting dynamical systems in their own right. In these 
devices, current can be used both to create and to de- 
tect vibrational excitations. Clear evidence of phonon 
excitations induced by single electron tunneling has been 
reported in a number of different systems, including 
semiconducting phonon cavities^ molecules^ and sus- 
pended carbon nanotubesiiSiii At finite bias electrons 
tend to drive phonons out of equilibrium; signatures of 
non equilibrium phonon distribution were observed in a 
suspended carbon nanotubeiiS 

On the theoretical side, NEMS are often described 
with simple phcnomenological models involving a sin- 
gle electron device coupled to an harmonic oscilla- 

£ 0r 12.13.14.15.16.17.18.19.20.21.22.23.24.25.26.27.28 :n — jthin 

these simple models, many peculiar features such as neg- 
ative differential conductancefi 2 ^*^ shuttling instabil- 
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and strong mechanical feedback^^have been 
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predicted in the case of an undcrdamped oscillator 
is then a physically relevant question whether the vi- 
brational energy is reduced by relaxation processes in- 
duced by coupling to an external environment^ or rather 
because tunneling itselfii2ii2iS£ Up to now, theoretical 
works have focused mostly on the case of negligible relax- 
ation, taking the opposite case of strong relaxation as a 
reference term. Significant differences between these two 
cases have been found both for weak and strong electron- 



and F was always found to be larger or equal 
However, Fano factors below this limit in the 



phonon (e-ph) couplingiSSiSi 

Many recent theoretical works have focused on the 
study of current noise on NEMSi 19 i 20 i 21 i 22 i 23 i 24 i 25 In par- 
ticular, the Fano factor F, which is the ratio between the 
zero frequency component of the noise and the average 
current, has proven to be very sensitive to the e-ph inter- 
action and to the details of the phonon distributioni^SiSi 
A giant enhancement of Fano factor (F ~ 10 2 — 10 3 ) 
has been predicted for strong coupling and negligible 
phonon relaxation^ In the opposite limit of strong relax- 
ation, i.e. when the phonons arc thermally distributed, a 
NEMS behaves essentially as a single electron transistor 
(SET)^i Shot noise in SET has been extensively stud 
iec 

to 1/2 

single electron tunneling regime have been predicted in 
more complicated systems. For istance, coupling to in- 
ternal degrees of freedom^ 3 , can induce F slightly below 
1/2 (F ~ 0.45). A strong suppression of the noise has 
been predicted for the quantum shuttleiS 5 ^ 2 ^ In this case, 
very low values of the Fano factors stem from an highly 
ordered charge transfer mechanism given by strong cor- 
relations between charge and mechanical motion. 

In this work we discuss how intermediate phonon 
relaxation influence the transport properties of a SET 
coupled to a mechanical oscillator. We focus on the 
sequential tunneling regime and we adopt a rate equation 
to describe the dynamics of the system. This approach 
is justified when the characteristic frequency of the os- 
cillator is much bigger than the tunneling ratefiSiiiSSiS 7 . 
which is the typical experimental situation. We find that 
finite relaxation rate affects the dynamics in a highly 
non trivial way. Both current and noise can be either 
enhanced or suppressed by relaxation, depending on the 
e-ph coupling and on the considered voltage range. In 
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particular, for voltages higher than a certain critical 
value, the Fano factor can be even suppressed below 1/2. 
This suppression is observed in a completely incoherent 
regime as a consequence of the interplay between vibra- 
tion assisted tunneling and direct relaxation of different 
vibrational states. 

The paper is organized as follows. The model Hamil- 
tonian is defined in Sec. [HI while in Scc. lllll wc introduce 
the rate equation and the formal expressions for the cur- 
rent and the noise. In Sec. IIVI numerical results for the 
current and Fano factor are presented: in particular, the 
suppression of the Fano factor is discussed in detail for 
a wide range of parameter. Finally, analytic expressions 
for the current and the Fano factor are derived within a 
toy model employing few phononic states. 



II. MODEL 

In several experimental realizations, either using litho- 
graphically defined quantum dots^£ molecules^ or nan- 
otubesjifliii electron transport is dominated by single 
electron tunneling— In this regime, the system is essen- 
tially a SET coupled to an harmonic oscillator. Describ- 
ing the SET as a single electronic level, the Hamiltonian 
of the system is H s = H n +H b + H njb where- 20 ! 21 (ft = 1) 

H„ =en, (1) 
H b =cu (tfb+ 1/2), (2) 
Hn, b = Aw (6t+6)n. (3) 

The operator n = cftd represents the occupation number 
of the single level, whose energy e = e(V g ) can be tuned 
with the aid of an external gate voltage V g . Vibrational 
excitations are created by and their ground state is 
defined as the zero-phonon state when n — 0. The fre- 
quency of the oscillator loq can range from the hundreds of 
MHz of a nanometrical cantilever— to a dozen of THz in 
the case of molecular devices or suspended nanotubesAifi 
The dimensionless parameter A in the coupling term H 7lib 
represents the strength of the e-ph interaction. For ex- 
ample, A ~ 1 was reported for the Ceo devices^ and for 
suspended carbon nanotubes jii while values of A between 
0.4 and 3 have been found in different C140 samples^ 

The SET is coupled to external leads by a tunneling 
Hamiltonian 

H t = ^( c L d + dtc fc,"), (4) 

where the operators c\ create electrons with momen- 
tum k in lead a = 1,2. The leads are described as non 
interacting Fermi liquids with 

Hleads = ^ £k,ac\ a Ck,a (5) 
k,a=l,2 



and their chemical potential can be shifted by a bias 
voltage V. For simplicity, in the following we will as- 
sume symmetric voltage drops and symmetric barriers 
t\ = t<i = to- 

Finally, the oscillator is coupled a dissipative environ- 
ment that we describe as a set of harmonic oscillators^ 

H env =5^w J -(o]o,- + l/2), (6) 

3 

H b ,env = ^2xjVj(<4 + «j)( &t + b )- ( 7 ) 
3 

Here are the creation operators of the bosonic bath 
modes. The environmental coupling is usefully charac- 
terized by its spectral function 

JH=2i^^%^ J ). (8) 

3 

III. RATE EQUATION 

The eigenstates of H s can be written as \n, I), where 
n denotes the occupation of the single level and I the 
phonon number. The coupling to the leads and to the 
environment induces an energy broadening of these eigen- 
states. If this broadening is the smallest energy scale of 
the problem, a perturbative treatment for H t and H b env 
is appropriate and a master equation for the reduced den- 
sity matrix of the system can be derived in the sequential 
tunneling regime— At lowest order, the reduced density 
matrix is diagonal in n but may still be off-diagonal in I 
because of the e-ph coupling H n b . 

In the following we will consider the case where loq 
is much larger than the bare tunneling rate = 
2-Kvt\ (with v the density of states of the leads). In 
this "diabatic" regime, which is the typical experimen- 
tal situation)2i2ii£ the elements of the density matrix 
that are non-diagonal in phonon number become neg- 
ligible jiSiiiSSiSl Then, the master equation reduces to a 
rate equation for the occupation probabilities P n i of the 
state \n, I) 

jl-^ni = ^ y^[-Pn'i'Iq P->j" ~ Piling "_>;" ] 

+ Y,\ P nl'T r itl-PnlT r i%]- (9) 

V 

The coefficients T a represent the tunneling rates 

through the a-th barrier while T^,, are the relaxation 
rates. 

In order to evaluate such rates, it is convenient to 
eliminate the coupling term H n b from H s by means of 
a canonical transformation. Due to the coupling term 
Hb,env, this transformation must include both the oper- 
ators of the oscillator and those of the environmen1ii£ 

O = e An Oe- An , A = - b) - 2 K ^x 3 {a] - a s ), 

3 
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where 



1 - 4 Eix|^/"o' 

The total Hamiltonian is transformed into 



(10) 



H — H n + Hb + Hi ea ds + H env + H t + Hb 



where 

Ht= Yl to(cl a e~ A d + Se A c k , a ). (11) 

fc,a=l,2 

and H n = en with s = e — As the energy of the 

SET is renormalized by a factor proportional to A 2 , this 
represents the relevant parameter for the e-ph interac- 
tion. 

The transition rates can now be easily calculated using 
Fermi golden rule, giving rise to tunneling rates propor- 
tional to ig, and relaxation rates, which depend on Xj- 

The relaxation rates represent transitions between vi- 
brational excitations without change of the electronic 
state ((3- 1 = k B T) 
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(12) 



where J(uiq) is the spectral density of the phonon bath 
Eq. (|SJ), evaluated at the frequency of the oscillator. 
Treating Hb_ env at second order allows only transitions 
between neighboring states (i.e. \V — l\ = 1). Transitions 
with \V — l\ > 1 can be included considering different 
relaxation mechanisms 

The charge transfer rates are induced by Ht ■ Assuming 
the electrons in the leads are at equilibrium with their 
chemical potential, one obtains the following expressions 



'-a l- 
L a I- 



rWxmUW-i)), (13) 
r^x n [i -f a (u; (i-i'))}, (14) 



where f a {x) = f(x + e — Su a ), f(x) is the Fermi func- 
tion and 8fi a = (— l) a+1 eV/2 is the shift of the chemical 
potential of the leads induced by the bias voltage. The 
coefficients Xw are given by 



\{n, l\e 
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(A 2 )| 



(15) 



where l K = min{/,/'}, £> = max{/,/'} and Lf(x) is a 
generalized Laguerrc polynomial. These terms are called 
Franck- Condon factors and are well known from molec- 
ular spectroscopy^ The effect of the e-ph interaction on 
transport is two-fold: on one hand it suppresses the ef- 
fective tunneling rate (because of the factor e _A ), on 
the other it induces a non-trivial dependence on the 
phononic indices 1,1' . Up to moderate e-ph coupling 
(A 2 < 1), transitions which conserve or change slightly 
I have the largest amplitude and those between states 



with low vibrational number are dominant. Vice versa, 
the latter are exponentially suppressed for A 2 3> 1, while 
transitions which change I considerably become favored. 

Within the rate equation approach, the current and 
noise can be evaluated by means of standard tech- 
niques^2i2i It is convenient to adopt a matrix formalism 



and write the rate equation as 

\P) = M\P), 



(16) 



where the vector \P) = {P n i\ represents the time depen- 
dent occupation probability distribution. Calling \p( st >) 
the stationary solution of Eq. I|ltj|) . the steady current 
through the a-th barrier is 



(/„> = e(l\l a \P^), 



1,2 



(17) 



where (1| = (1,1,..., 1), e is the the charge of the 
electron and I a is a matrix representing all the possi- 
ble transitions through the considered barrier (X a )V"'y = 

(-l) Q+1 (?i - n') r «™'^z™- Following Ref.21, the zero fre- 
quency component of the noise is given by 

S a A°) = -^e 2 {l\{5X a M- l 51 l3 + 81pM- l 5X a )\P {st) ) 
+2 e 2 ^(l||Z Q ||P^), (18) 

where ST a = T a — (I a )/e. Because of charge conserva- 
tion, the steady current and the zero frequency current 
correlators are independent of barrier index: {I a } = /, 
Sa,p(0) = S. In the following we will mostly refer to the 
Fano factor F = S/2eI. 



IV. RESULTS 

A. Full solution of the rate equation 

The dynamics of the system is characterized by two 
competing time scales: the average time spent by an elec- 
tron in the dot r e ; and the phonon relaxation time T p h- 
If T e i 3> T ph , the vibrational excitations tend to relax be- 
tween each tunneling event to the thermal Bosc distribu- 



tion P, 



(eq) 



(1 



). In this limit, charge and 



(eq) 



vibrational degrees of freedom decouple P n i = P n P\ 
and the dynamics of the system reduces to an effective 
two-state sequential tunneling process^ The analytic 
expressions for current and noise are well known^S and, 
for ksT -C eV, are respectively given by 



j(eq) 



e — 



Aeq) 



f 2 



fi 



(Ti 



(19) 



Here f 1 = r<°> £i ai/iM and f 2 = T^^a^l - 
/ 2 (— Iojq)] are the renormalized rates for tunneling in and 
out of the dot and a; arc Poissonian weight factors^ 
ai = 6(l)e~ x X 2l /l\. In this case the smallest possible 
value of the Fano factor is F^ = 1/2. 
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Vice versa, if r e ; <c r p h the tunneling electrons drive 
the vibrations out of equilibrium and peculiar features 
such as negative differential conductance (NDC)i^i^i 
and super-Poissonian shot-noise2i have been predicted. 

In our model, a rough estimate of t^ 1 is given by t p 7 1 = 

r(o) e - 



i.e. by the effective transparency of the barrier 

ph 



set by the e-ph coupling, while r ^ is determined by 



environment spectral density r~ h = J[ujq). It is useful 
to define a dimensionless parameter for the relaxation 
strength 



w = JM/r^ 



(20) 



In terms of w, the condition for equilibrated phonons 
T e i ^ T p h reads w 3> cxp(— A 2 ). It is then evident that 
the e-ph coupling defines a characteristic scale for relax- 
ation: the stronger is the coupling, the more sensitive is 
the system to phonon relaxation. 

This is reflected by the stationary phonon distribu- 

For increasing relaxation 



tion P, 



(st) 



Pi 



(st) 



01 



Pi 



(st) 



1 / 



strength w, P^' tends monotonically to P^ eq ^ but with 
A 2 -depending speed (see Fig^). For strong e-ph cou- 
pling, the phonon distribution is narrow already in the 
non-relaxed case 2 ^ w = and it reaches equilibrium for 
values of w which are sensibly smaller than for weak A 2 . 




6 0,60; 



FIG. 1: Stationary phonon probability distribution P^ s *' for 
different values of w. Upper panels: A 2 = 0.4, lower: A 2 = 7. 
The rightmost panels (w = oo) represent the thermal Bose 
distribution P ; (e<7) = e~ mu °(l - eT^ w °). Other parameters: 
V = 3 loo, e = and k B T = 0.02 uj . 



Since P£ st ' converges monotonically to pj- eq ^ for grow- 
ing w, one expects most of the features of the non equi- 
librated case to be washed out by increasing relaxation. 
This is particularly evident in the case of the giant Fano 
factor observed at low voltages for strong interaction 
(A 2 ^> I), which is strongly suppressed even by weak 
relaxation (see Fig. [2J . This behavior can be easily un- 
derstood observing that F>1 depends dramatically on 



the non equilibrium distribution of the vibrational exci- 
tations induced by tunnclingi 21 i 28 In fact, for large A 2 
transitions between low lying phonon states are expo- 
nentially suppressed (see Eq. I|15[l h Therefore, the main 
contribution to the current comes from high excited vi- 
brational states (states with large I) but at low voltages 
the occupation probability of those states is strongly sup- 
pressed ^i22i These conditions leads to avalanches of tun- 
neling processes which, in turn, are responsible for the 
huge values of Fi 2 i2& Direct phonon relaxation inhibits 
this mechanism reducing even further the occupation of 
states with large I and, consequently, both the current 
and the Fano factor are strongly suppressed. For very 
strong relaxation (w — > oo), F — > 1/2 as one would ex- 
pect for equilibrated phonon on resonance (e = 0). 
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FIG. 2: Fano factor as a function of voltage for A 2 = 16 and 
for different values of the relaxation strength w; red w = 0, 
green w = 0.1, blue w = 1, magenta w = 10. Dashed line, 
F — 1/2. Other parameters: e = and ksT — 0.02o;o. 

Similarly, relaxation has a destructive effect on NDC 
(not shown) as this is also a consequence of the peculiar- 
ity of the nonequilibrium phonon distribution induced by 
tunneling itselfiia2L 

One could be tempted to conclude that considering ex- 
plicitly the effects of relaxation simply results in an "in- 
terpolating" behavior between the opposite limits of no 
relaxation and thermally distributed phonons. However, 
we find that finite relaxation rate can induce unexpected 
features. 

Let's first consider the case of moderate coupling A 2 = 
3. In Fig. we plot the Fano factor function of 
voltage for different values of w and for e = 0. It appears 
that F has a non systematic dependence on id: it can be 
cither enhanced or suppressed by relaxation depending 
on the considered voltage range. For eV < 6luq it is 
always F > 1/2. In particular, for eV < 2wo it is F = 1/2 
as the tunneling electrons cannot excite vibrations and 
the system behaves as an ordinary single level. More 
interestingly, for eV > Qluq, relaxation can suppress F 
even below 1/2. 

It is worthwhile to stress that we are dealing with a 
single electronic level in the sequential tunneling regime, 
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FIG. 3: Fano factor as a function of voltage at A =3 and 
for different values of the relaxation strength w; red w = 0, 
green w = 5, blue w = 12, magenta w = 100. Dashed line, 
F = 1/2. Other parameters: e = and fcsT = 0.02ojo- 

coupled to a harmonic oscillator with a completely inco- 
herent dynamics. Therefore, this unexpected suppression 
of the Fano below 1/2 can only be ascribed to the in- 
terplay between vibration-assisted tunneling and direct 
relaxation of the phononic excitations. Indeed, relax- 
ation induces a tendency to ordered transfer of electrons 
thought the SET via emission-absorption of phonons. 




FIG. 4: Density plot of the Fano factor as a function of bias 
V and relaxation w for different values of A. In all the panels: 
dark gray F < 1/2, medium gray F — 1/2 (indicated by 
the arrow in the color map ) and light gray F > 1/2. The 
white line, corresponding to F = 1/2, represents wt(V). The 
black line in the 4th panel delimits the region where noise 
is superpoissonian, F > 1. Other parameters: e = and 
k B T = 0.02^0- 

From a numerical analysis, it emerges that this pecu- 
liar behavior can be found for any value of the e-ph cou- 
pling. In particular, we observed that it exists a voltage 



threshold set by the e-ph interaction 

eVt(A) = 2lu int[A 2 ] (21) 

such that, for V > Vt(A), relaxation larger than a cer- 
tain threshold value wt (V, A) suppresses the Fano factor 
below 1/2. This is shown in Fig. which represents 
a grayscale plot of the Fano factor in the (V, ui)-plane, 
for different values of A 2 . The white contour line corre- 
sponds to F — 1/2 and separates two different regions 
in the (V, u;)-plane: the one to the right of the contour, 
where F < 1/2 and the other one where F > 1/2. In 
other words, the white line denotes Wt as a function of 
V at given A 2 . The threshold voltage V t (A) corresponds 
to the position of the vertical asymptote of wt(V, A). For 
A 2 < 2 the critical voltage coincides with the onset of 
vibration assisted tunneling eVt(X) = 2loq; vice versa for 
strong e-ph coupling (A 2 3> 1), V*(A) becomes very large 
and this is why F is always higher than 1/2 in Fig. [21 

The minimal value assumed by the Fano factor f m j„ 
depends itself on the e-ph coupling (see Fig. |SJ). For 
weak coupling, F m i n differs only slightly from 1/2. For 
stronger coupling (A 2 > 1) it decreases logarithmically 
and it only reaches the value F m i n ~ 0.4 for consider- 
ably strong interactions. Note that each point in Fig. [5] 
corresponds to different values of voltage and relaxation 
strength, as the position of F rn i n in the (V, w)-plane de- 
pends on A 2 . The inset shows the voltage V m in where 
the minimum is found. 
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FIG. 5: Main panel: F m in as a function of A 2 . Each point 
corresponds to different values of w and V. Inset: voltage 
Vmin where the minimum is found as a function of A 2 . 

Finally, let's observe that for A 2 > 2, the threshold 
voltage 14(A) corresponds to the onset of the transition 
I ; — * int[A 2 ]. In this case V*(A) is a characteristic volt- 
age also for the current which can be either suppressed 
or enhanced by relaxation depending on V being smaller 
or larger than Vt(A) (see Fig.|SJ. Relaxation contributes 
to populate the low lying phonon states and then, at 
low voltages, it inhibits the current as the transitions be- 
tween those states have exponentially suppressed rates. 
However, for V > Vt(X) the transition I : — > int[A 2 ] 
is allowed and, as it correspond the greatest Franck- 
Condon factor^ it gives a substantial contribution to 
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the current. In this case relaxation has the opposite ef- 
fect and it sustains the current "feeding" the population 
of the vibrational ground state. For V ~ Vt(\), these 
two mechanisms coexist and, consequently, the current 
depends only weakly on relaxation (see inset in Fig. [BJ) . 

This observation fits nicely what is reported in litera- 
turei^SiSi In fact, for A 2 < 2 the critical voltage is smaller 
than the energy required to have phonon-assisted tunnel- 
ing and then the current is enhanced by phonon relax- 
ation at any voltage, consistently to what was observed 
in Ref^i Vice versa for very strong e-ph coupling the en- 
hancement of the current due to relaxation can be hardly 
seenSi as Vt(A) shifts to very large voltages. 




eV/uJ 

FIG. 6: Current as a function of voltage for A 2 = 3 and for 
different values of the relaxation strength w; red w = 0, green 
w = 1, blue w = 10, magenta w = 100, cyan w = oo. Other 
parameters e = and fcsT = 0.02ojo- Inset: zoom of the 
plateau around eV = 5lJq. Current in units el^ '. 



B. Toy model 

To get some insight in these results, we focus on the 
low voltage region eV < 4u>q and we consider a toy model 
with only four accessible state, i.e. n = 0, 1 and I = 0, 1. 
In this case analytical expressions for current and noise 
can be derived. For sake of simplicity, we report only the 
solutions on resonance (e = 0) and at zero temperature 



I = e r<°> 

and 

F= l - 6 -l 



x m , a , T , „ X 01 {w + 2X ol - 

—+e(eV-2uj ) 

2 2(w + 2X i) 



A) 



(22) 



eV - 2oj )X iA[w 2 + w(2X 01 - A) - X 01 A] 



(w + 2X i) 2 K 
[w(X 01 



X n 



(23) 



A) 



where A = Xqq — Xn and K 
Xoi(2X i+2X 11 +A)]. 

From Eq. I|22|) it is easy to show that the current is 
an increasing function of w only for A > (that is, for 
A 2 < 2, see Eq. |JTSJ|). Vice versa, for A < (A 2 > 2) the 



current decreases for increasing relaxation, in agreement 
with what was previously observed. Moreover Eq. I|23|) 
tells that A > (A 2 < 2) is the necessary condition to 
have F < 1/2 in the region 2wo < eV < 4^0- In fact only 
in this case, it exist a threshold value for relaxation 



2w t (X) = A-2X Q1 + JA 2 



(24) 



such that for w > w t the Fano factor is smaller than 
1/2. For stronger e-ph coupling A 2 > 2 (A < 0), it is 
always F > 1/2. This confirms the numerical estimate 
eVt(A) = 2wo as the threshold voltage for any A 2 < 2. 

Despite the coarseness of the model, Eq. (|23f) ac- 
cords qualitatively with the exact numerical solution for 
eV < 4u>o (see Fig.[7J). The agreement is reasonably good 
even in the case of weak e-ph coupling, where the phonon 
distribution is mostly broadened^ and one expects the 
four state approximation to be more inaccurate. A bet- 
ter agreement can be obtained considering a six state 
model with n — 0, 1 and I = 0, 1, 2 but, in this case, 
the analytic solutions become quite cumbersome and we 
don't report them here for simplicity. The agreement 
of the four states model with numerical result suggests 
that F < 1/2 rather depends on the interplay between 
relaxation and vibration-assisted tunneling, than on the 
possibility to access an high number of vibrational states. 
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FIG. 7: Fano factor as a function of relaxation strength w for 
eV — 3u>o and for different values of the e-ph interaction. Red 
curve: exact numerical result; blue curve: result for the four 
states model - Eq. (ll'Mt : green curve: result for a six states 
model with n = {0, 1} and I = {0, 1,2}. Other parameters 
e = and k B T = 0.02^ - 



V. CONCLUSIONS 

In conclusion, we have investigated the effects of 
direct phonon relaxation on the shot noise of a SET 
coupled to a mechanical oscillator. For increasing 
relaxation strength, the occupation probability distribu- 
tion of the states of the system evolves monotonically 
towards thermal equilibrium. In contrast, we found a 
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non-monotonous behavior of the Fano factor, which can 
be suppressed even below 1/2. This relaxation-induced 
tendency to order of the electronic transfer through 
the dot is unexpected, since we are dealing with an 
oscillator with incoherent dynamics coupled to a SET 
in the sequential tunneling regime. The onset of this 
behavior is discussed as a function of relaxation, e-ph 
interactions, and external voltages. We have found that 
for any value of the e-ph coupling, a critical voltage 
Vj(A) exists such that for V > V t (X) a suppression of 
the Fano factor below 1/2 is possible. At low voltages, 



these results are qualitatively predicted by a four states 
toy model. In this work we have focused on the case 
of a symmetric device on resonance, which is the most 
favorable to analyze the suppression of the Fano factor. 
However, qualitatively analogous results can be obtained 
for asymmetric barriers and (or) e =/= 0. 
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